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RADIAL REFLECTION DIFFRACTION TOMOGRAPHY ( 
[0001] The United States Government has rights in this invention pursuant to 
Contract No. W-7405-ENG-48 between the United States Department of Energy and 
the University of California for the operation of Lawrence Livermore National 
Laboratory. 

RELATED APPLICATION 
[0002] This application claims the benefit of U.S. Provisional Application No. 
60/474,861, filed May 30, 2003, entitled, "Radial Reflection Diffraction Tomography/' 
which is incorporated herein by this reference. 

BACKGROUND OF THE INVENTION 

Field of the Invention 

[0003] The present invention relates in general to imaging, and more particularly 
to an imaging method and apparatus employing Radial Reflection Diffraction 
Tomography. 
Description of Related Art 

[0004] Intravascular ultrasound (IVUS) imaging provides a method for imaging 
the interior of blood vessel walls. In standard acoustical techniques, a catheter with a 
rotating ultrasound transducer is inserted into a blood vessel. The transducer 
launches a pulse and collects the reflected signals from the surrounding tissue. 
Conventional ultrasonic imaging systems use B-mode tomography or B-scans, 
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wherein images are formed from the envelope of the received display echoes 
returning to an ultrasonic transducer as brightness levels proportional to the echo 
amplitude and by assuming straight ray theory (i.e., geometrical optics). The 
brightness levels can then be used to create cross-sectional images of the object in the 
plane perpendicular to the transducer image. However, such images typically suffer 
from the consequences of ray theory of sound propagation, which does not model its 
wave nature. 

[0005] A circumferential scan can be made by either rotating a single transducer 
(mechanical beam steering) or by phasing an array of transducers around a 
circumference (electronic beam forming). Typically, one ultrasound pulse is 
transmitted and all echoes from the surface to the deepest range are recorded before 
the ultrasound beam moves on to the next scan line position where pulse 
transmission and echo recording are repeated. When utilizing B-scan, the vertical 
position, which provides depth of each bright dot is determined by the time delay 
from pulse transmission to return of the corresponding echo, and the horizontal 
position by the location of the receiving transducer element. 
[0006] Although B-scan IVUS images can be utilized to detect plaque and 
characterize its volume, the classification of plaque types by ultrasound is very 
difficult. Conventional B-scan images utilizes scattering, which, in turn, depends on 
the acoustic impedance dissimilarity of tissue types and structures. Although hard 
calcifications in some plaque can be detected using such a mismatch, the similarity 
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in the acoustic properties of fibrous plaque and lipid pools prevents direct 
identification. Consequently, the size of the fibrous cap cannot be accurately 
estimated. 

[0007] Diffraction tomography has additionally been applied to medical imaging 
problems for a number of years, usually in a standard circumferential through 
transmission mode. Furthermore, improved vascular images have been provided by 
utilizing time domain diffraction tomography, a technique capable of accounting for 
the wave propagation of the transmitted acoustic waves in addition to redundant 
information from multiple angular views of the objects imaged. A related B-mode 
approach that incorporates spatial compounding has also been employed to provide 
improved vascular images through multiple look angles. 

[0008] Background information on rotational IVUS systems are described, for 
example, in U.S. Pat. No. 6,221,015 to Yock. Background information on phased- 
array IVUS systems are described, for example, in U.S. Pat. No. 6,283,920 to Eberle et 
al., as well as U.S. Pat. No. 6, 283,921 to Nix et al. Multi-functional devices have been 
proposed in other areas of vascular intervention. For example, U.S. Pat. No. 
5,906,580 to Kline-Schoder et al., describes an ultrasound transducer array that may 
transmit signals at multiple frequencies and may be used for both ultrasound 
imaging and ultrasound therapy. Therapeutic ultrasound catheters, are described, 
for example, in U. S. Pat. No. 5,725,494 to Brisken et al. and U.S. Pat. No. 5,581, 144 to 
Corl et al., which describes another ultrasound transducer array that is capable of 
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operating at multiple frequencies. However, none of the above devices and 
associated techniques from the above cited patents, are suited for rapid 
identification of objects, such as, but not limited to, vulnerable plaque or objects 
recessed in a bore hole, in accordance with the principles of the present invention. 

SUMMARY OF THE INVENTION 
[0009] The present invention is directed to a wave-based imaging method, which 
includes: directing predetermined energy waves radially outward from within an 
interspace and receiving scattered energy waves from one or more objects. The received 
data are processed to produce images of the objects, wherein the processing includes 
application of a wave-based algorithm that can map an angular location and a plurality 
of frequency parameters of the received scattered energy waves to construct images of 
the one or more objects. 

[00010] Another aspect of the present invention is directed to a wave-based imaging 
method that can be utilized to characterize a plaque, which includes: inserting a 
catheter having a longitudinal axis and a distal end into an artery, wherein the catheter 
further includes a single transmitter disposed about the distal end of the catheter and a 
receiver aperture having a plurality of receivers additionally disposed about the distal 
end of the catheter, wherein the transmitter and the receiver aperture is capable of 
rotating up to 360 degrees about the longitudinal axis of the catheter. As part of the 
method, one or more predetermined energy waves are directed radially outward from 
the single transmitter and radial scattered energy waves are received in a 
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predetermined imaging mode by the receiver aperture. The received scattered energy 
waves results in collected data capable of being processed to produce images of plaque 
from the surrounding artery walls, wherein the processing includes application of a 
wave-based algorithm that can map an angular location and a plurality of frequency 
parameters of the received scattered energy waves to construct the images and 
determine the risk of rupture and/ or thrombosis. 

[00011] Another aspect of the present invention is directed to a wave-based imaging 
method that can be utilized to characterize a plaque, which includes: inserting a 
catheter into an artery, directing one or more predetermined energy waves radially 
outward and receiving one or more radial scattered energy waves from a distal end of 
the catheter; collecting a radial scattered tomographic data baseline of the artery's 
tissue; measuring an applied external pressure to the artery; obtaining a deformation 
radial scattered tomographic data set of the artery's tissue after application of the 
external pressure; and processing the radial scattered tomographic data baseline and 
the deformation radial scattered tomographic data set to produce a final image 
indicating elasticity of the artery to characterize the imaged plaque, wherein the 
processing includes application of a wave-based algorithm that can map an angular 
location and a plurality of frequency parameters of the received scattered energy waves. 
[00012] A further aspect of the present invention is directed to a wave-based 
imaging apparatus, which includes a flexible substrate having a longitudinal axis 
and a distal region and one or more elements disposed on the distal region and 



5 



Attorney Docket No. IL-10883 

Express Mail Certificate No. EV283924528US 



capable of directing one or more predetermined energy waves radially outward and 
receiving one or more radial scattered energy waves from one or more objects. The 
received scattered energy waves are capable of producing images of one or more 
objects by processing a collected data set, wherein the processing includes 
application of a wave-based algorithm that can map an angular location and a 
plurality of frequency parameters of the received scattered energy waves. 
[00013] A final aspect of the present invention is directed to a wave-based 
imaging apparatus that includes a Hilbert space inverse wave (HSIW) algorithm 
that can map an angular location and a plurality of frequency parameters of said 
received reflected diffracted energy waves so as to characterize plaque within a 
living vessel. 

[00014] Accordingly, the present system and method employs desired Radial 
Reflection Diffraction Tomographic techniques to determine the state and location of 
buried wastes, to track plumes of underground contaminants of materials, to 
determine the state of materials residing in waste drum barrels or weapons, to 
evaluate nondestructively parts having existing access holes (e.g., automobile parts), 
and for identifying potentially life threatening vulnerable plaque buildup on living 
vessel walls. 

BRIEF DESCRIPTION OF THE DRAWINGS 
Fig. la illustrates a basic multimonostatic mode configuration that 
includes a single transducer rotating about a fixed center. 
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Fig. lb illustrates a basic multistatic mode configuration that includes 
a fixed annular array of outwardly directed transducers. 

Fig. lc illustrates a basic multistatic mode configuration that includes a 
rotating aperture. 

Fig. 2 shows a conventional IVUS catheter. 

Fig. 3a shows a conventional IVUS catheter inserted into a diseased 

artery. 

Fig. 3b illustrates the RRDT geometry of the present invention when a 
catheter is inserted into a diseased artery. 

Fig. 4 illustrates RRDT non-destructive evaluation within a bore hole. 
DETAILED DESCRIPTION OF THE INVENTION 
[00015] Referring now to the following detailed information, and to incorporated 
materials; a detailed description of the invention, including specific embodiments, is 
presented. 

[00016] Unless otherwise indicated, numbers expressing quantities of ingredients, 
constituents, reaction conditions and so forth used in the specification and claims are 
to be understood as being modified by the term "about." Accordingly, unless 
indicated to the contrary, the numerical parameters set forth in the specification and 
attached claims are approximations that may vary depending upon the desired 
properties sought to be obtained by the subject matter presented herein. At the very 
least, and not as an attempt to limit the application of the doctrine of equivalents to 
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the scope of the claims, each numerical parameter should at least be construed in 
light of the number of reported significant digits and by applying ordinary rounding 
techniques. Notwithstanding that the numerical ranges and parameters setting forth 
the broad scope of the subject matter presented herein are approximations, the 
numerical values set forth in the specific examples are reported as precisely as 
possible. Any numerical value, however, inherently contains certain errors 
necessarily resulting from the standard deviation found in their respective testing 
measurements. 
General Description 

[00017] The present invention employs inverse wave techniques to reconstruct 
images of a medium surrounding a physical probe in a plane perpendicular to an 
axis of rotation in a radial reflection configuration, i.e., in a multimonostatic or a 
multistatic arrangement disclosed infra, wherein one or more transmitting and 
receiving elements, more often at least about 15 of such elements, e.g., 
transducer(s), are at a fixed radius and designed to collect scattered fields, e.g., 
reflected and diffracted fields. Such a radial reflection diffraction tomography 
(RRDT) technique is based upon a linearized scattering model to form images 
given the disclosed physical transmitter and receiver configurations and the 
mathematical method, i.e., a Hilbert-based wave algorithm, utilized to invert the 
scattering collected fields. 
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[00018] As example embodiments, the multimonstatic and multistatic probes 
can be mounted at the end of a flexible substrate, such as a catheter or snaking 
tube that can be inserted into a part with an existing access hole or a medium 
(e.g., an artery) with the purpose of forming images of the plane perpendicular to 
the axis of rotation. By applying the Hilbert space inverse wave (HSIW) 
algorithm of the present invention to the collected data of such multimonostatic 
and multistatic probes, radial reflected diffraction tomographic images are 
readily obtained. 
Specific Description 

[00019] Fig, 1(a) shows a basic multimonostatic conceptual arrangement of the 
present invention, wherein a single energy source element 1, such as a transducer, 
can operate as both source and receiver (as denoted by T/R, to indicate transmitter 
and receiver) at multiple spatial locations. At each angular location along the 
illustrated dashed circumference, as denoted by the directional arrow, energy source 
element 1 can launch a primary field wave (not shown) and receive a reflected 
scattered field wave (not shown). Such an arrangement often requires a spectrally 
wide band frequency diverse source capable of producing frequencies from about 1 
kHz to about 3 THz (Electromagnetic frequencies), often from about 100 Hz to about 
10 GHz (Acoustic frequencies), more often between about 20 MHz and about 60 
MHz (Acoustic frequencies), to provide spatial diversity so as to form images of a 
surrounding medium. 



9 



Attorney Docket No. IL-10883 

Express Mail Certificate No. EV283924528US 



[00020] Fig. 1(b) shows an example conceptual multistatic mode embodiment, 
wherein a plurality of fixed energy source elements 2, e.g., transducers, are arranged 
as an annular array, generally designated by the reference numeral 20. In succession, 
along for example, the illustrated directional arrow, each energy source element (for 
example, the element denoted by the letter T to indicate a transmitter) is capable of 
launching a primary field wave (not shown) and a backscattered field wave (not 
shown) is measured on all the remaining elements (as denoted by the letter R, 
indicating the remaining fixed elements are operating as receivers). 
[00021] Fig. 1(c) illustrates a beneficial multistatic conceptual arrangement that 
includes a plurality of energy source elements 4 configured in a rotating sub- 
aperture 6, as denoted by the bi-directional arrow, formed by a single transmitter, as 
denoted by the letter T to indicate a transmitter, and surrounded by multiple 
receivers, as denoted by the letter R. At each angular location, as denoted by the 
single directional arrow along the illustrated dashed circumference, transmitter T 
can launch a primary field (not shown) and a backscattered field (not shown) is 
measured on all receivers R. 

[00022] When operating in a reflection mode as disclosed herein, the mathematics 
applied to the collected data operate beneficially to image objects because the range 
resolution of the reconstructed image is proportional to the number of frequencies 
used in the reconstruction. Under the Hilbert space inverse wave algorithm, 
increasing the number of frequencies and transducers, increases the complexity of 
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the reconstruction, the size of the intermediary data files, reconstruction time, and 
computer memory requirement. Thus, the trade-off between computer resources 
and resolution is a consideration. Nonetheless, the techniques employed in the 
present invention are beneficial even at acoustic frequencies from as low as about 
100 Hz to as high as about 10 GHz. Such lower frequencies allow the disclosed 
embodiments to additionally be employed in borehole type of applications, such as, 
but not limited to, characterizing underground contamination plumes or waste in 
contamination barrels. 

[00023] For either the multimonostatic or multistatic example embodiments, the 
planar reconstruction of the imaged object(s) requires that the one or more collected 
measurements map a pair of spatial variables (i.e., angular location and incident 
source frequency) of a physical object into the angular location and frequency 
parameters of the measured field. 

[00024] An exemplary application of the present invention is in the 
characterization of vulnerable atherosclerotic plaque. Arthrosclerosis is a condition 
where the arteries are obstructed by the buildup of deposits, "intravascular plaque" 
(IVP), on the inside of arterial walls and such a buildup of deposits can lead to what 
is known to those of ordinary skill in the art as cerebrovascular disease, which is the 
third leading cause of death and the leading cause of major disability among adults. 
Plaque grows as a fibrous tissue encapsulating a lipid pool and as the plaque grow it 
may incorporate calcium. Of particular concern is vulnerable or unstable plaque 
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because of the possibility of such plaque becoming inflamed and unexpectedly 
rupturing. Stable or non-vulnerable plaque, typically includes a thick layer of 
fibrous tissue of about 800 microns but is not life threatening and can be treated 
slowly. A thin fibrous cap of typically up to about 300 microns covering a pool of a 
soft lipid core typically characterizes vulnerable plaque. When such a cap is 
disrupted, the thin cap is compromised and the lipid deposited into the artery can 
produce adverse effects, such as thrombus formation, strokes and death. 

[00025] Fig. 2 shows a conventional catheter, generally designated as reference 
numeral 200, for intravascular tissue characterization, such as atherosclerotic 
vulnerable plaque. Such a catheter 200, typically has an elongated flexible substrate 
202 with an axially extending lumen 204 through which a guide wire 206, and /or 
various other devices or other instruments can be delivered to a region of interest. 
An ultrasonic transducer assembly 208 is provided at the distal end 210 of the 
catheter, with a connector 214 located at the proximal end of the catheter for 
transducer manipulation and processing received transducer signals. Transducer 
assembly 208 can comprise a plurality of transducer elements 216 arranged in a 
cylindrical array centered about a longitudinal axis 218 of the catheter for 
transmitting and receiving ultrasonic energy. Adhesive (not shown) and or an end- 
cap (not shown) can be applied to transducer assembly 208, and lumen 204 to 
protect such elements from the surrounding environment. Transducer assembly's 
208 individual elements (not shown) and conductive acoustical backing (not shown) 
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are often mounted on the inner wall of elongated flexible body 202 operating as the 
flexible substrate. 

[00026] Fig. 3a illustrates a typical IVUS method using such a catheter 200, as 
shown in Fig. 2a. In such a conventional application, catheter 200, having a 
transducer assembly 208 that can launch an energy wave as a primary field (as 
denoted by the letter F) is inserted typically non-centered into a nominally circular 
diseased artery 302. Around a wall 304 of artery 302 is a fibrous collagen plaque 306. 
A lipid pool 308 can reside inside such a fibrous structure, wherein when a fibrous 
cap 310 of plaque 306 separating lipid pool 302 from the blood (not shown) within 
artery 302 is more than about 800 /xm thick, plaque 306 is characterized as stable. 
However, in cases where cap 310 is less than about 300 /im thick, such a plaque is 
characterized as vulnerable, and is more likely to rupture and /or thrombosis. 

[00027] The present invention utilizes the disclosed RRDT approach for improved 
intravascular applications such as characterizing plaque as discussed above, and 
incorporates various aspects of the method of utilizing a probe, such as, but not 
limited to, the catheter as shown in Fig. 3a. However, such catheters 200 and similar 
probes known to those of ordinary skill in the art typically show angular overlap for 
beam processing, which results in loss of valuable image information of one or more 
objects of interest within a surrounding medium. The present invention overcomes 
such processing by incorporating novel improvements of the transmitters and 
receivers, by utilizing frequencies between about 20 MHz (Acoustic) and about 60 
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MHz (Acoustic), and by utilizing RRDT techniques of the present invention as 
discussed herein. Such novel embodiments accounts for phase, amplitude, and 
beam diffraction, to recover not only such loss of valuable image information 
information but to further enhance the imaging capabilities of the invention by 
providing images with improved lateral resolution of the acoustic absorption and 
sound speed. 

[00028] Fig. 3b shows the geometry incorporated by the RRDT method of the 
present invention. Fig. 3b shows a cross-sectional view of a catheter 200, having an 
outer diameter between about 0.25 mm and about 5 mm, being inserted into an 
artery 302, having a surrounding plaque 306 that includes a cap 310 and a lipid pool 
308. Inserted into artery 302 is a non-centered catheter 200, which includes a 
transducer assembly (not shown) that can be disposed about the distal end of 
catheter 200, as disclosed in the present invention, with a radial location specified by 
r 0 = /? o (cos0 o ,sin 0 O ) , where R 0 is the catheter 200 probe radius, a constant. At each 
angular location, 9 0 , transducer assembly 208, as shown in Fig. 3a, launches a 
primary field F radially outward (as denoted by the letter r) into a medium, such as 
the blood (not shown) and surrounding tissue in this example, and the transducer 
arrangement, as disclosed in the present invention, can measure a reflected scattered 
field (not shown) having, for example, at least up to about 90 degrees of angular 
content from one or more objects, such as the linings of cap 310 that overlies lipid 
pool 308. 
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[00029] As another example embodiment, the RRDT method and apparatus of the 
present invention can be combined with elastography to gain further insight into a 
surrounding medium's elastic properties and provide further information in the 
determination of characterizing plaque as vulnerable or stable. Generally, the 
contrast in elastic properties between a lipid pool and a fibrous cap is evident. By 
utilizing elastography, the elastic properties of a vessel wall can be obtained by 
observing a deformation of the vessel due to an external pressure, such as the 
pressure produced by a heart. Such a change in the arterial pressure due to the 
pumping action of the heart produces a measurable deformation of the tissue 
surrounding the vessel. Such a deformation can be measured by tracking a motion 
of patterns in successive intravascular scans as disclosed by the present invention. 
By knowing the arterial pressure and the measured deformation, the present 
invention can recover elastic properties of the surrounding tissue. From such elastic 
properties, one can further characterize the surrounding tissue to predict plaque 
composition. 

[00030] Fig. 4 illustrates a further beneficial embodiment, wherein the present 
invention can be utilized for non-destructive characterization (i.e., RRDT imaging) in 
applications other than for intravascular RRDT imaging. As shown by the example 
cross-sectional underground view of a borehole 404 in Fig. 4, a flexible substrate 400 
or snake-like tube having a transducer assembly 402 similarly configured like the 
intravascular RRDT application discussed above, can be lowered into bore hole 404 
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so as to image a site using RRDT techniques. Such an arrangement can launch a 
primary field (denoted by the letter F) and receive diffracted energy waves having 
frequencies often between about 100 Hz and about 300 Hz, to determine the state of 
buried wastes, such as waste within a radioactive waste drum barrel 410 or a 
biohazardous container, and /or to track a plume 412 of underground contaminants. 
In a similar manner, disclosed probes herein, can be inserted into waste drum 
barrels 410, or weapons (not shown) or any part having an existing access hole, such 
as, but not limited to, an automobile engine, and determine the state of the part or 
material. 

Hilbert Space Wave Inversion 

[00031] Hilbert spaces are spaces constructed using vectors. Specifically they 
define vector spaces where sets of vectors in the space "add up" to another 
vector, an analog to Euclidean space where measurements can be added to result 
in another valid measurement. Hilbert spaces are particularly useful when 
studying the Fourier expansion of a particular function. In the Fourier transform, 
a complex function describing a waveform is re-expressed (transformed) into the 
sum of many simpler wave functions. A Hilbert space describes the "universe of 
possible solutions" given one particular such function. The Hilbert space inverse 
wave (HSIW) algorithm of the present invention enables an inverse for any 
multistage or multimonostatic geometry with any combination of sources, 
receivers, and frequencies. 
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[00032] In a radial reflection device of the present invention, such as an 
intravascular ultrasound probe having an outer diameter between about 0.25 mm 
and about 5 mm, or a probe configured to non-destructively characterize buried 
wastes (e.g., tracking plumes of underground contaminants of materials), evaluating 
the state of materials residing in waste drum barrels or weapons, or to non- 
destructively evaluate parts with existing access holes (e.g., automobile parts), 
acquired data are collected at discrete angular locations. Such angular locations are 
denoted by: 

(1) i^s^Ccos^.sin^) 

for transmitter locations, where 0 n -nA0 src ior n = 0,l -,N src -1, where 
N src 27tl A0 , and A0 src is the source angular increment. 

Similarly, receiver locations are given by: 

(2) R r m ^R o (cos0 m ,sm0J 

where 0 m =mA0 rcv for m = 0,l—,tf JfC -l, where N rcv 2nl A0 rcv , and A0 rcv is the 
receiver angular increment. 

For each source, configured receiver(s) can record a backscattered field as 
a time series that can be digitized for processing. Discrete Fourier transforming 
the time series data result in the spectrum of one or more measured wave forms 
at discrete frequencies. The forward scattering equation under the Born 
approximation with both spatial and frequency diversity is given by: 
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(3) iC'(/M>^/) = ^^ 



where 0) n l = 0,1, N f -1 are the discrete frequencies and is the number 

of frequencies in the pulse band width. 

The HSIW as disclosed'herein interprets Equation (3) as a mapping 
from a continuous object space to a discrete measurement space. The object 
space is the physical (x,y)space of the object function. The measurement space 
includes discrete angles and temporal frequencies at which the scattered data 
are collected. The scattering operator projects the object onto the 
measurement space. The forward propagation or projection kernel is defined 
as: 



where H(r) is a J = (N src x N rcv x N f ) element column vector, and P(a> { ) is the 
incident pulse spectrum. Mathematically, the projection is represented as an 
inner product between the object function and the kernel via: 



where D is a J element column vector, and where each element represents a 



(4) 



n\r)^P{0) l )k 2 o {(o l )G{R r m ,r,(o l )G{r,R , n ^ l ) t 



(5) 



D = \drU\r)o(r)^ai,o), 



particular source, receiver, and frequency combination. Symbolically, the 



forward scattering operator, K, is denned as: 



(6) 



K[]^idrU\r)[]. 
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The HSIW method of the present invention is employed to derive an 
inverse of the operator as shown in equation (6). The singular value 
decomposition (SVD) of K is given as: 
(7) K = USV\ 

where the columns of U form an orthonormal set of column vectors, w y , 

which span a measured data space, and the components of V form an 
orthonormal set of vectors, v y (r) , which span an object space. S is a diagonal 

matrix of singular values, o . . It is emphasized that the u j are complex 

column vectors where as the v y (r) are complex functions of r. The set of 

normal equations for such a singular system are: 



(8) Kv j (r) = a j u j , 

(9) *\.=<7 y v.(r), 

(10) KK'u j =a j Kv j (r) = <j>u j/ 

(11) K'Kvj(r) = (TjK'ujir) = a> ; (r) , 



The inversion method of the present invention estimates the object 
function of equation (5) given measured data in D. Such an inversion 
incorporates expanding the object function in terms of v ; (r) : 

(12) 5(r) = l>;v.(r), 

7=0 

where the a. are constant coefficients to be determined. Substituting the 
object expansion into equation (5) results in: 
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7-1 7-1 



(13) D = \dr n*(r)5> y v,.(r) = \drYV(r) Vj (r) , 

Applying the definition of the K operator in equation (6) to equation 
(8) yields an expression for the integral of equation (13), 

(14) Kv j =\drYT(r)v j (r) = a j u j/ 
which reduces equation (13) to: 

y-i 

(15) D = J5>/r y u,, 



Using the orthogdnality of the u. vectors, the unknown a j is solved as 
follows: 



y-i y-i 



(16) u]D = Y. a F^ u i = Y* a i a A = <*Pi > 

7=0 7=0 

resulting in: 

(17) a,.&. 



The adjoint of the forward scattering operator, K and the singular 
values and singular vectors, a u jf and v ; (r) are now required. First, the 
following inner product equation defines the adjoint, 
(18) < M ,tfv> = <tf + M ,v>, 
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Using the definition of the forward scattering operator from equation 
(16) results in: 

(19) w + IdrYT(rMr) = Yl\r))v{r), 

By comparing the right hand sides of equations (18) and (19), the 
following definition of the adjoint of the forward scattering operator is 
obtained: 

(20) *t[.] s [.]. n ». 

The a . and u } are determined by solving the eigenvalue equation of 

equation (10) formed by the outer product of the forward scattering operator 
with its adjoint. Explicitly, the outer product is represented by: 

(21) (ldrir(r)U T (r))u j =^u j/ 

which is a / x / eigenvalue equation of the form Ax = Ax . The ri(r) vectors are 
known analytically and can be evaluated numerically. It follows that the 
elements of the outer product matrix can be computed numerically and the 
resulting system solved numerically for the a) and u j . Given these and using 
equation (19) to solve for v ; .(r) results in: 

(22) v(r)= ±n T (r)u j . 

°> 

Substituting equations (17) and (22) into equation (12) yields the final 
expression for the reconstruction: 
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As described above, the Yl(r) vectors of equation (4), and outer 
products and eigenvalues of equation (21) are computed numerically. The 
measurement system of the analytically described invention only measures 
part of the scattered field due to the aperture and the loss of the evanescent 
field information and accordingly, some of the eigenvalues, cr y 2 , are close to 

zero. Those eigenvalues and their corresponding eigenvectors determine the 
rank of the outer product matrix, and they must not be used in the 
reconstruction of equation (23). Thus, in the method of the present invention, 
a Best Rank N approximation is used to select the number of singular 
values/vectors. A ratio is computed as follows: 

N-\ 

(24) R(N)=iS— , 

7=0 

where the singular values are assumedly arranged from smallest to largest: 
al < of < cj)_ x . Plotting R(N ), the point at which the function starts to rise 
rapidly is graphically identified. The index of the singular value at which this 
occurs is labeled as J 0 . With this value determined, a final reconstruction is 
as follows: 
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(25) d(r)=^-^YV(r)u j u)D 

i=h °) 

[00033] The HSIW as disclosed herein is flexible in that it allows any transducer 
configurations of the present invention and any number of frequencies to be used in 
forming such a final reconstruction. 

[00034] Changes and modifications in the specifically described embodiments can 
be carried out without departing from the scope of the invention, which is intended 
to be limited only by the scope of the claims. 
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